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- - . It is shown that recent criticism by C. R. Hagen questioning the validity of stress tensor treatments 

' of the Casimir energy for space divided into two parts by a spherical boundary is without foundation. 
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I. INTRODUCTION 



< 

oo 

In a recent paper, Hagen criticizes conventional, well-established techniques for computing the Casimir stress on 
spherical, or in general, curved boundaries In our view this paper is full of misconceptions and incorrect assertions, 
■ and sheds no light on the subject. The point seems to be that the Green's function method for calculating the Casimir 
stress on a spherical shell, as originally proposed in Ref. ||^, is incorrect, yet, mysteriously, manages to arrive at the 
correct answer The simpler, and correct conclusion, would be, on the contrary, that that method is valid. As we 
' will demonstrate, such is the case. 
: Three errors are claimed to have been made in such calculations: (i) the discontinuity of the stress tensor Tij has 
] been used to calculate the force per unit area on a spherical surface; (ii) poles are neglected in the rotation to Euclidean 
space; and (iii) the incorrect outgoing wave boundary conditions have been imposed on the Green's functions. On 
the contrary, all of the procedures used in the criticized papers are, in fact, correct. 

Let us first address the point (iii). The causal or Feynman propagator is used both in the Milton et al. papers [§J^-j|] 
and in Hagen's paper. Hagen seems not to appreciate the simple identity obeyed by the free photon propagator: 



S ■ D^{v t) = r —e-^-' f — - r dk sinfcre-^^l*! - H ^e"''^* — e^l^l'' (1) 



where the second form is obtained by first integrating over to, the third by integrating over k. The last form explicitly 
. , shows that the casual propagator corresponds to outgoing-wave boundary conditions.P] A pedagogical derivation of 
Eq. (|l|) and other properties of the causal propagator is given in Problem 31.9 in Ref. p^. 
" " ' As to point (i): There can be no doubt that the discontinuity of T„ across the spherical surface represents the 
stress per unit area exerted on that surface. As shown in the Appendix of Milton and Ng |15| ], this is equivalent to 
the form given by the stationary principle of electrostatics, the force density being f — —^E'^e for a dielectric body. 
(See for example, Ref. [0.) Moreover, in most of Milton's papers, consistency was checked explicitly by rederiving 
the formulae by at least two different methods, either from the stress on the surface, or from the energy obtained by 
integrating Tqq over all space. In fact in the original paper |^ a third independent variational method was used. The 
stress tensor formulation of the problem is further validated in Sec. 3. The problem that Hagen addresses seems to 
be the result of using undefined formal expressions [for example, the sum in his Eq. (8) does not exist] so statements 
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^It is also possible to use the retarded Green's function, which follows from a statistical mechanics argument In general 

the connection between the causal and the retarded Green's dyadics is Q'r+(r, r',a;) = sgn(a;)SjG^(r, r', tj). The same results 
are obtained. 
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about coefficients vanishing in the limit of the enclosing sphere becoming infinite are not meaningful, for one must 
regulate the expressions before a limit can be taken. 

The final objection (ii) concerns the existence of poles in the complex frequency plane. Indeed such poles occur 
in the lower half frequency plane for the exterior contribution, when the exterior volume is infinite (no large exterior 
sphere). But exactly this issue was dealt with in the original paper on p. 395, and will be explained in detail in 
the next section. 

So none of Hagen's objections appear to be valid. Of course, it is not incorrect to enclose the sphere in question by 
a much larger sphere to keep the eigenvalues real, but it is not necessary, and the simplified procedure is correct. 

There are many other minor misstatements and incomplete remarks that render Hagen's paper objectionable. The 
recent experiments |p^ , while admittedly based on a somewhat different geometry than parallel plates, indeed confirm 
the theory at the 5-1% level because the correction to a spherical lens is easily done. Moreover, there can be little 
room for scepticism about the reality of Casimir forces, in view of the confirmation of the closely related Lifshitz 
theory [p^j9|JlO[| , and the demonstration of the equivalence of van der Waals forces with the Casimir force between and 
within dilute bodies . The Casimir force for cylindrical and spherical geometries has recently been confirmed by 
several authors using completely different zeta function techniques |2^. These, and other methods, also give finite 
results for other geometries besides the "electromagnetic sphere," namely for cylinders, for scalar and fermionic modes, 
for dielectric bodies when e/i = const, etc. We could go on, but our point is taken. 



II. ON THE LEGITIMACY OF ROTATING THE INTEGRATION CONTOUR 



Here we readdress in more detail the issue of "rotating the contour of frequency integration." As we will see, what 
is involved is a bit more sophisticated: It is a Euclidean transformation. The argument as given in |^ applied to the 
Green's functions, either inside or outside the sphere of radius a, for which either r,r' < a, or r, r' > a. In fact, as 
employed in that reference and elsewhere, the "rotation" is actually applied to the expression for the Casimir energy, 
or the stress at the surface. To perform the transformation, it is essential that both inside and outside contributions 
be included. We can, however, consider the scalar (or TE) modes separately. Those give for the force per unit area 
on the sphere (where we have not included a large external sphere — its inclusion does not modify the argument) 



^i+l /2 i^^) '^1 + 1/2 (^") 



(2) 



where k = (For the derivation of Eq. (||) see [^.) Here the sign of the additive constant has been reversed, in 
effect, by adding a contact term, a term proportional to <5(t), so that the resulting frequency integral be convergent. 
Now it may be verified that the integrand in Eq. (|^) has the following analytic properties in the complex variable 
C = ka: 

• The singularities lie in the lower half plane or on the real axis. Consequently, the integration contour C in a; 
lies just above the real axis for a; > 0, and just below the real axis for u; < 0. 

• For 3C > 0, the integrand goes to zero as 1/|CP- (This is a weaker condition than specified in ||^.) This 
convergent behavior is the result of including both interior and exterior contributions. 

Then we may write the stress on the sphere as 

/ = / 5(|c^l), (3) 

Jc 2vr 



where the integrand satisfies the dispersion relation 

TT* J-oo C 



gi\L0\) = - / 2 ■ 9iC), (4) 



because the singularities of g(C) occur only for 3^ < 0. Now we can carry out the to integral in Eq. (g) to obtain 

1 „-i|ClkL 



Finally, we re-write the result in Euclidean space by making the Euclidean transformation i\T\ —> |t4| > 0, so that we 
have the representation 



2 



2ici j^^2^ki+e- 

Thus the EucHdean transform of the stress is 

f^J'^-'l ^^'"'^'ami). (7) 

In effect, then, the Euchdean transformation is given by the recipe uj iki, |w| — > i\ki\, r — > ir^. In particular, the 
force per unit area given for a massless scalar field (0) is transformed into the expression 



JE = 



1=0 



,.,^(2^ + 1)^/ ^ye-W.^^+.;^ + l , (8) 



where we have adopted dimensionless variables, y = ua, 6 = r^/a 0, x = \y\. It is now straightforward to evaluate 
this expression. The result is as (surprisingly) first given in 1994 / = 0.0028168/(47ra"'). Thus a careful analysis 
shows that the methodology in Ref. ||] is correct, and that, in fact, as long recognized, a finite result is only achievable 
if the difference between Trr{r = a — e) and Trr{r = a + e) is taken. 



III. ON THE LOCAL INTERPRETATION OF THE STRESS TENSOR FOR A CURVED BOUNDARY 



A. Theory 

Let us readdress the issue about the relation between the change in the stress tensor T"^^ across an arbitrary 
surface, and the resulting force density /" on the surface. In Hagen's paper it is argued that the well-known relation 
between /" and a change in T„ across a planar surface in Cartesian coordinates is lost when one is dealing with 
curved surfaces, and apparently even when one is using non-Cartesian coordinates independent of the presence or 
non-presence of curved surfaces. Hagen bases his reasoning on the assumption that Vq.T"'' = holds everywhere. 
However, it is an almost trivial fact that this is not the case when sources for the electromagnetic field are present. 
It is in particular straightforward to show that one in general has 

v„r"^ = -F^^ Ja - -Z*^ (9) 

where is the source current, F°'^ is the electromagnetic field strength tensor, and is the four- vector force density 
acting on the material sources in the system. Since this relation is covariant we conclude that the interpretation of 
/" as a force density holds in any situation. This interpretation is in particular independent of the explicit system of 
coordinates employed or whether one is dealing with a curved spacetime. One can always locally recover the usual 
Minkowski spacetime results by transforming to a local vierbein. 

The total force on the source is simply the volume integral of the force density in some three-dimensional 
spacelike hyper-surface S. The total force vector field on a simply connected volume V in S with boundary dV is 
most simply computed using Stokes' theorem 

= ( = - / VaT"'5 = - / '^"^"'^ • (10) 

Jv Jv JdV 

In this expression represents a unit spacelike outward normal to V. We will now specialize to the case when we 
are dealing with a thin conducting spherical shell, and when Maxwell's vacuum field equations hold everywhere away 
from the shell. Thus V represents the shell, and dV its inner and outer surfaces. Clearly, in general the total stress F 
in the direction of a unit vector field u" is = Jy f^up. Due to time- independence and the exact spherical symmetry 
of the shell, we choose it" to point in the radial direction, and we find that the total outward stress on the shell per 
unit solid angle is 

J^^ = r2T'^'-(r_) -r^r"(r+). (11) 

In this expression the components of the stress tensor are to be evaluated at the inner and the outer radii of the shell, 
which are denoted by r_ and r+ respectively. In the limit when the shell is infinitely thin (r_|_ —^ r_) with radius a 
we thus find that the force per unit area is 
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(12) 

The same result can also be obtained by integrating the divergence of the stress tensor directly, using the fact that 
Tr(TQ/3) — 0. This should prove once and for all that Hagen's critical discussion about the use of the change in Trr 
as a direct measure of the stress on a spherical shell is, from the theoretical side, completely wrong. Let us now turn 
to the experimental side. 



B. Experiments 

We wish to present two counterexamples from optics which demonstrate that Hagen's claim about the use of T„ in 
calculating the force on a spherical shell is incorrect in this regard: The predictions obtained from local stress tensor 
calculations are in accordance with experimental evidence. 

1. The first example concerns the oscillations of a small water droplet illuminated by a laser pulse. In 1988 
Zhang and Chang |^l| took photographs of water droplets (50 /Ltm radius) at various instants after illumination by 
a A=0.60 /im dye laser pulse of 0.4 /xs duration. In the case of a 100 mJ pulse, the droplet began oscillating with 
surface elevations of up to 30 per cent of the equilibrium radius. Since a spherical drop acts as a lens, the maximum 
elevations were observed, as expected, at the rear section. 

From a theoretical point of view this problem is a hybrid problem: the natural way to proceed is to calculate the 



electromagnetic surface force density using Eq. (12) as if the surface of the sphere were at rest. Thereafter one can 
calculate the hydrodynamic response of the sphere by means of the Navier-Stokes equation. The problem was studied 
by Lai tt al. [E2j , and by Brevik and Kluge p3[ , for various polarizations in the linear approximation for the surface 
elevation. Wave optics was required. Although the linear approximation may appear rough at the rear section, it is 
clear from the figures that this kind of theory is in reasonable agreement with the observations. This example, as far 
as we are aware, is one of the very few cases where the observations of the radiation pressure are local. The viability 
of the expression ([l^) is quite evident. 

2. Our second example is about the optical forces observed on microparticles in an evanescent laser field. There is 
a striking experimental demonstration of this effect, by Kawata and Sugiura [ p4| : The particles are lifted up from the 
dielectric interface and moved parallel to the interface. Although this effect is not a local effect such as above, one has 
nevertheless, in order to describe it theoretically, to make use of the same formula (O) as before. The theory of the 
effect has been worked out by Almaas and Brevik |2^, Lester and Nieto-Vesperinas |p6||, and others. The agreement 
between wave theory and observations is reasonable, again justifying the legitimacy of Eq. (E3). 
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